A multi-wavelength analysis of Spitzer selected Coma Cluster 
galaxies: star formation rates and masses 
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ABSTRACT 

We present a thorough study of the specific star formation rates for MIPS 24/im selected galax- 
ies in the Coma cluster. We build galaxy spectral energy distributions using optical (u',g',r',i',z'), 
Near-infrared (J,H,K S ), an d Mid to Far-infrared (IRAC and MIPS) photometry. New and archival 
spectra confirm 210 cluster members. Subsequently, the total infrared luminosity, galaxy stellar 
mass, and specific star formation rate for the members are determined by measuring best-fit 
templates. Using an array of complementary diagnostics, we search for contaminating AGN, but 
find few. We compare obscured star formation rates to unobscured rates derived from extinction- 
corrected Ha emission line measurements. The agreement between these two values leads us to 
conclude that there is no evidence for an additionally obscured component. In our spectroscopic 
sample, complete to 80% for r'< 19.5, we find that all starbursts are blue and are dwarfs, having 
masses <10 9 M Q . Examining the location of these starbursts within the cluster, we confirm that 
there is a lower fraction in the cluster core. 



Subject headings: galaxies: clusters: individual (Coma Cluster) 
iesdundamental parameters (luminosities, masses) 



infrared: galaxies - galax- 



1. Introduction 

Galaxy clusters are dynamic systems in which 
tens to hundreds of galaxies are plowing through 
hot intracluster gas. In the core of the clus- 
ter, the galaxy density is the highest, and star 
forming galaxies preferentially avoid these regions 



evolutionary history of galaxies at the outskirts 
compared to those at the cluster core. However, 
the study was based on a spectroscopic sample 
and less deep than other studies. More recent 
studies, however, have revealed a segregation of 
galaxy activity with local density: Galaxies at the 



( Kennicutt 


1983 


Dressier et al. 


1999 


Gavazzi 


2009 


). However, at the outskirts of the cluster 



cluster center are oldest (Poggianti et al. 2001 



many galaxies show starburst activity detected at 



Mid-Infrar ed (IR) ( |Koyama et al.||2010| |Gallazzi 
et al.|2009 ) and radio ( Heckman et al.|1985 ) wave- 
lengths, as well as in optical emission lines (Rines 
eTatpiO^ . 

Coma, at z^0.023, is our nearest rich galaxy 
cluster. As such, it allows for a level of detailed 
study currently unavailable for more distant clus- 
ters. Mobasher et al. (2003) studied the R-band 



Smith et al. 


2008 


), with the fraction of blue, spi- 


ral, (Terlevich et al.|2001 


Aguerri et al.|2004 Bi- 


viano et al. 


2002 


, and dwarf (Odell et al. 


2002 


Gavazzi et al. 2010) galaxies increasing with dis- 



tance from the core. In the radio, a luminosity 
function has been well studied down to a very low 
surface brightnesses (Miller et al.| 2009 1. These 



luminosity function, finding no difference in the 



faint radio galaxies are vast in number, have high 
specific star formation rates (sSFRs), strong opti- 
cal lines, and appear to have undergone a recent 
burst of star formation, prior to their quenching. 
Metallicity gradients in cluster galaxies have also 



1 



been noted, and linked to the age segregation, as 
well as to pressure confinement from the deep po- 



2001 



tential at the cluster core (Poggianti et al 
Carter et aLfl2002| ). 

A scenario in which the galaxy activity is 
quenched during infall toward the dense cluster 



2009) 



core is supported by galaxy alignments (Adami 
et al.|200"9 1 , alpha element gradients ( Smith et al 



2004). 



sizes of spiral disks (Aguerri et al 
and the correlation of post starburst galaxy po- 
sitions with the location of X-ray substructures 



(Poggianti et al. 


2004) 


hajan et al. 


20111. 1 



population of recently 
quenched elliptical galaxies would also explain the 
observations of |Clemens et al.| ( |2009[ ) who studied 
the color-magnitude diagram of Coma at Mid-IR 
wavelengths finding that 68% of ellipticals were 
truly passive, but 32% showed evidence for very 
recent activity. 

This picture is consistent with results from 
cosmological simulations that show the build-up 
of massive clusters through accretion of infalling 
galaxies and galaxy groups. The numerical simu- 
lations of Bekki ( 1999) show higher star formation 
rates (SFRs) in the accreted galaxies, in agreement 
with recent observational results. Gall azzi et al.l 



(20091 found that most of the obscured star for- 
mation in the complex Abell 901-Abell 902 comes 
from galaxies between and at the edges of the two 



clusters. Our work on Abell 1763 ( jEdwards et al. 



2010||Fadda et al.||2008| ) revealed a -12 Mpc fila- 
ment of Mid-IR emitting galaxies connecting this 
cluster to the nearby Abell 1770. For that case, 
the fraction of star forming filament galaxies was 
found to be double that in the cluster regions. 



Work by Jenkins et al. (2007), Bai et al. (2006) 



and Mahajan et al. (2010) has shown that the 



dwarf star forming galaxy population in Coma is 
quenched in the high density core. However, the 
dwarf galaxies do not necessarily have very intense 
starbursts, therefore, a study of the specific star 
formation rates is important. This requires an ac- 
curate method for determining galaxy stellar mass 
determinations. Thus, of particular importance to 
this study are our new Near-IR observations as 
they bracket the stellar mass bump, allowing for 
accurate stellar mass measurements. In this pa- 
per, we combine these original Near-IR photome- 
try with public Mid- Infrared photometry from the 
IRAC and MIPS instruments aboard the Spitzer 



Space Telescope, in addition to optical photom- 
etry from the Sloan Digital Sky Survey (SDSS). 
We gather optical spectra from the SDSS as well as 
from our own campaign. With these data we build 
spectral energy distributions (SEDs) of Coma clus- 
ter members. Star formation rates are calculated 
from the total infrared luminosity derived from 
best-fit template galaxy SEDs. 

In the next section, we present the observations 
and reduction methods for the deep Near-IR J, H, 
and K s data which was collected from the Palomar 
WIRC instrument, as well as the Spitzer MIPS 24, 
70, and 160/im and IRAC 3.6, 4.5, 5.8, and 8.0^m 
data which we have re-reduced. We also include 
the optical photometry from the SDSS Data Re- 
lease 7, the available spectra from that program, 
as well as 338 spectra of 24/im selected galaxies 
we gathered from the HYDRA instrument on the 
WIYN telescope. In Section 3, we present the 
catalog of 24/im selected galaxies in Coma, along 
with the associated optical and infrared fluxes and 
rcdshifts. We release our fully reduced images 
and spectra to the public via the NASA/IPAC In- 
frared Science Archive (IRSA) and the NASA Ex- 
tragalactic Database (NED). These data allow us 
to determine the amount of AGN contamination, 
measure stellar mass, and derive the total infrared 
luminosities required to calculate specific star for- 
mation rates (SFRs), which we present Section 4. 
We compare these obscured SFRs to unobscured 
rates measured from nebular emission lines in Sec- 
tion 5, where we also compare other proxies for 
the obscured SFR and discuss the specific SFR as 
a function of location. We conclude, in Section 
6, that the starburst galaxies are blue dwarfs and 
that there is a paucity of these in the dense cluster 
core. This is in line with the hypothesis that star- 
burst galaxies have been recently accreted onto the 
cluster and will be quenched as they make their 
way to the core. 

2. Observations and Data Reduction 

In this section, we describe the observations and 
data reduction. Table [I] lists our observational 
dataset and Figure [T] shows the coverage of the in- 
frared data and the location of spectroscopic mem- 
bers. 

The photometric data was acquired from re- 
cent WIRC Near-IR and archival Mid and Far- 
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Table 1 
Observations 



INST 


A cc „ t (/im) 


ID 


Date 


Timc(min) 


Coverage(' 2 ) 


Depth 


FWHM (") 


IRAC 


3.6, 4.5, 5.8, 8.0 


r3859712 


2004 Jul 04 


160.1 


3250 


1.5, 4.0, 3.9, 3.8 /iJy 5ct 


1.7,1.7,1.9,2.0 


IRAC 


3.6, 4.5, 5.8, 8.0 


r3859456 


2003 dec 18 


110.9 








IRAC 


3.6, 4.5, 5.8, 8.0 


r3859968 


2004 Jul 04 


160.1 








MIPS 


24, 70, 160 


r4740096 


2004 Jun 22 


169.6 


11070, 10530 


0.08, 32.0, 130 mjy 5er 


5.9, 16, 36 


MIPS 


24, 70, 160 


r4740352 


2004 Jun 22 


169.6 








MIPS 


24, 70, 160 


rl8315776 


2007 Jan 19 


170.3 








MIPS 


24, 70, 160 


rl8316032 


2007 Jan 19 


170.3 








MIPS 


24, 70, 160 


r-18315520 


2007 Jan 19 


170.3 








MIPS 


24, 70, 160 


1-18316288 


2007 Jan 19 


170.3 








WIRCJ 


1.250 




2008 Apr 23 


250.0 


8480 


19.5 mag 


1.2 


WIRCH 


1.635 




2008 Apr 25 


187.0 


8480 


18.5 mag 


0.9 


WIRCK S 


2.150 




2008 Apr 24 


247.0 


8480 


17.0 mag 


1.3 



IR Spitzer imaging (Pis: Fazio, Rieke, Kitayama), 
optical catalogs from the SDSS DR7, and radio 
source catalogs from the Faint Images of the Ra- 
dio Sky at Twenty-Centimeters (FIRST). Because 
SED fitting is more accurate when the MIPS 24/im 
datapoint exists, we include only galaxies which 
are detected at 24/tm to at least the 3c level. The 
MIPS data field of view covers the full core of the 
cluster, (0.2 - 0.5)xr 2 oo (where r 2 oo is the radius 
within which the cluster density is 200 times the 
mean density required to close the universe). This 
same region has been mapped in J, H, and K s us- 
ing the WIRC instrument on the Palomar 200 in 
telescope. 

We also take advantage of spectra from SDSS 
DR7 which we add to our own data from WIYN. 

2.1. WIRC Near-IR Observations 

The Near-IR magnitudes are driven by old 
stars, which comprise the bulk of the galaxy mass 
and are therefore paramount in our later deriva- 
tion of the galaxy stellar mass. 

The Near-IR photometry was gathered at the 
Palomar 200in telescope using the Widefield In- 
fraRed Camera (WIRC) during three nights of ob- 
servations between April 23-25, 2008. As detailed 
in Table 1, the nights were clear with seeing be- 
tween 0.9-1.3" and most of the region covered by 
Spitzer was imaged, as shown in Figure [l] Fields 
in J were observed in sequences of 40s, and for 30s 
in H. In K s , 3 coadds of 13s each were taken to 
avoid saturation. The total integration times for 
J, H, and K s are 250.0, 187.0, and 247.0 minutes, 
respectively. Dark frames were also obtained for 



each integration time. 

The same basic reduction technique was applied 
to all three bands using the pipeline of Tom Jar- 
ret Q A median dark frame was subtracted from 
the data frames, and correction terms for the flux 
non-linearity are applied. A median sky made of 
a maximum of 10 frames was calculated and sub- 
tracted from the data frame before flat fielding. 
By accounting for the rotational offset of the tele- 
scope using a list of known stars from the Two Mi- 
cron All Sky Survey (Skrutskie et al. 2006; here- 
after, 2MASS) the astrometry was measured and 
verified by eye. The background level and flux bias 
corrections were then subtracted from the frames. 



SWarp (Bertin & Tissier 2007) was used to mo- 



saic the frames together. A final flux calibration 
on the mosaic, relative to the 2MASS catalog, was 
performed. 

The Near-IR photometry reaches a depth of 
J~19.5, H~18.5, and K s ~ 17.0 with uncertain- 
ties of 1.9, 2.7, and 4.6%, respectively. Figure [2] 
demonstrates how this is roughly two magnitudes 
deeper than those available with 2MASS. 



SExtractor (Bertin & Arnouts 1996) was used 
to identify sources and measure their magni- 
tudes. To ensure we extracted both small and 
extended objects, we used a det-minarea of 5 and 
det_minthresh of 2.0. To ensure we found the 
smaller objects that often surround larger objects 
(ex. near the brightest cluster galaxy) we used a 
deblend_mincont of 0.00005. 



1 http: / /spider. ipac.caltech.edu/staff/jarrett/wirc/jarrett. html 
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MIPS 




13.07 13.0 12.93 

RA(J2000) 



Fig. 1. — Coverage of IR and optical data. The coverage of the infrared photometry. MIPS 24^m 
coverage is represented by the large dotted rectangle, IRAC 3.6 and 5.8/zm fields are shown as smaller dashed 
squares in black and grey, respectively. The black polygon shows the Near-IR coverage. SDSS galaxies with 
r' < 20.5 and MIPS emission are shown as points, and those that are spectroscopically confirmed cluster 
members are large black dots. The filled grey circle has a radius of 0.3 degrees and highlights the cluster 
core (0.2xr 2 oo)- The MIPS data cover 8480 arcmin 2 . Deeper 24/xm observations exist over the area covered 
by the grey rectangle. 
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14 16 18 
Magnitude 



Fig. 2. — Depth of Near-IR data. Our new Near-IR images reach a depth of approximately two magni- 
tudes deeper than those available from 2MASS. 
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2.2. Spitzer Observations 

We have collected archival photometry from the 
Spitzer Space Telescope (Pis: Fazio, Rieke, Ki- 
tayama). The data were taken in 2003, 2004, 
and 2007 and all sets were observed for at least 
1.8 hours (see Table [I). The MIPS instrument 
provides Far-Infrared images at 24/im, 70/xm and 
160/im. IRAC data cover the Mid-Infrared wave- 
lengths of 3.6, 4.5, 5.8, and 8.0/im, in Figure [T] 
they are illustrated as dashed lines and available 
only for the cluster core. 

Each photometric band is analyzed on its own 
to produce a catalog of galaxies. Dust enshrouded 
bright stars re-emit their energy at longer wave- 
lengths, and the 24/im flux approaches this peak 
of the spectral energy distribution (SED), allowing 
for an estimate of star formation rates. Although 
MIPS 70/im and 160/im observations go even fur- 
ther toward the SED peak, these detectors are 
not as sensitive as those at 24/im. Therefore, we 
choose to base our catalogs on the 24/im sources. 
Sources from the 24/im catalog are then matched 
to candidates from each of the other photometric 



bands (see Section 3.1 for details). The main steps 
in data processing are outlined below. More de- 



tailed descriptions can be found in Edwards et al. 



2010 


) and 


Fadda et al. 


(2006) 



2.2.1. MIPS 24nm Data Reduction 

The Basic Calibrated Data (BCDs) and calibra- 
tion files are downloaded from Leopard, pipeline 
version S17.0.4, which uses a factor of 0.045 to con- 
vert instrumental units to MJy sr _1 . These BCDs 
have already been corrected for dark subtraction, 
muxbleed, detector linearization, flat fielding, and 
cosmic rays. We redid the pipeline jailbar correc- 
tion after masking bright sources, applied a droop 
and superflat correction for all BCDs in each AOR, 
and apply a flat for each position of the cryo- 
genic scan mirror. We removed the zodiacal light 
and corrected the astrometry to match the SDSS 
source positions. 

The BCDs were mosaiced using MOPEX 



(Makovoz & Khan (2005); Makovoz & Marleau 



(2005)) and extended sources were extracted us- 
ing SExtractor. We ran SExtractor a second time 
on the residual image, after multiplying by 141.08 
to convert MJy sr _1 to /iJy (as listed in the MIPS 



Data Handbook]^] p. 31). We computed aperture 
fluxes within a radius of 10.3" as well as Petrosian 
magnitudes with 4 radii and a minimum radius of 
5.5". We applied a factor of 1.167 for the aper- 
ture correction [[] and divided by 0.961 to account 
for the color correction as listed in the MIPS Data 
Handbook. The data were combined from two dif- 
ferent programs, leaving the central region much 
deeper. Although we included this deep strip in 
the reduced images, we note that the bulk of our 
analysis is based on a spectroscopic sample which 
corresponds only to the depth of the full MIPS 
image. By measuring the noise in the mosaiced 
frames, we arrived at a 5a depth of 80/iJy in the 
24/im mosaic. The final catalog includes sources 
with a signal to noise ratio (SNR) >3 sources, 
computed using an optimal aperture radius. 

There is one source, NGC 4921, which is ex- 
tremely extended and has many star forming 
clumps at 24/im. SExtractor failed for this source, 
therefore, we calculated the photometry by hand, 
fitting a Petrosian radius to the portion of the 
light profile which encompasses the same region 
as seen at r'. 

2.2.2. MIPS 70fim and 160^m Data Reduction 

For the longer wavelength MIPS observations, 
we downloaded the raw frames (DCEs) from the 
Spitzer archive and use the GeRT0pipel inc which 
finds and calibrates the slope of consecutive data 
reads. The conversion factor for instrumental 
units into MJy sr _1 is 702 for the 70/im data, and 
47.4, for the 160/xm data. For the 70/im data, we 
used a second pass filtering technique to remove 
streaking in the final image. We identified the 
causal bright sources in a first pass, and repeat 
the filtering on the masked image. We applied an 
additional correction to the BCDs to remove latent 
stim artifacts by modelling and then removing the 
effect of the stimflash on the median of the BCDs 
that follow. 

Virtually all of the extragalactic sources at 
these longer wavelengths are point sources. We 
mosaiced the frames and conducted the extraction 
of the sources using APEX with a bottom outlier 
threshold of 2 and a top threshold of 3. For the 

2 http: / /ssc. spitzer. caltech.edu/mips/dh/ 
3 http: / /ssc. spitzer. caltech.edu/mips/apercorr/ 
4 http: / /ssc. spitzer. caltech.edu/mips/gert/ 
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70/zm data, we used an aperture radius of 15", 
a color correction of 0.918, and an aperture cor- 
rection of 2.244. For 160/im, we use an aperture 
radius of 20", a color correction of 0.959, and an 
aperture correction of 3.124. By measuring the 
noise in the mosaiced frames, we arrived at a 5a 
depth of 32.0 mJy in the 70/mi mosaic and 130 mJy 
in the 160/xm mosaic. 

2. 2. 3. IRA C Data Reduction 

For all four IRAC bands, we included correc- 
tions for the column pulldown, jailbars, zodiacal 
light, subtraction of bright stars and background 
level, as well as a treatment for the background 
droop near bright stars. The only alteration to the 



process discussed in detail in Edwards et al. (2010) 



is that for IRAC 5.8 and 8.0/j,m, we included a 
correction for the superflat before subtracting the 
bright stars, this removes the gradient inside each 
BCD. 

The IRAC data reach a uniform 5<r depth of 

I. 5, 4.0, 3.9, and 3.8/LtJy in channels 1-4, respec- 
tively. As with the MIPS data, there is a deeper 
patch of 137arcmin 2 in each channel. Although 
the area of the deeper observations is consistent in 
all channels, the center of the deeper data is not 
consistent. Nonetheless, our spectroscopic photo- 
metric cut prevents the use of the deeper central 
data strip. 

The IRAC sources were extracted using SEx- 
tractor and sources with close pairs (within a 
projected separation <7") were flagged. For 
IRAC 3.6/mi, 4.5/im, IRAC 5.8/un, and IRAC 8.0/xm, 

II, 10, 7, and 9%, respectively, of the sources are 
close pairs. If an IRAC source was within 1" of 
a star as determined by the SDSS star/galaxy 
separator flag, it was flagged as a star. For IRAC 
3.6/zm, 4.5/im, IRAC 5.8^m, and IRAC 8.0/im, 6, 
7, 7, and 6% of the sources were flagged as stars, 
respectively. 

We calculated aperture fluxes using diameters 
of 4,6, and 12". The aperture corrections we 
used follow those of the SWIRE technical note 



(Surace et al. 2005) which lists corrections for 
radii of 1.9, 2.9 and 5.8". Following their work, 
we divided the aperture fluxes in each band by 
the corrections at each aperture radius. The 
corrections for each aperture radius are 0.736, 
0.87, and 0.96, for IRAC 3.6/xm; 0.716, 0.87, and 



0.95, for IRAC 4.5//m; 0.606, 0.80, and 0.94, for 
IRAC 5.8/im; and 0.543, 0.700, and 0.940 for 
IRAC 8.0/j,m. To estimate the total flux we also 
quote Petrosian fluxes. These are based on aper- 
tures whose radius is the average of the azimuthal 
light profile out to the radius where the local sur- 
face brightness is a factor of 0.2 times the mean 



surface brightness within the radius (Petrosian 



1976). We used two Petrosian radii and a mini- 



mum radius of 3", and applied the aperture cor- 
rection for extended sources as explained on the 
IRAC calibration webpage[^ We quote the uncer- 
tainty of the fluxes as those given by SExtractor, 
known to be underestimates as they only repre- 



sent the error from Poissonian noise (Valentinuzzi 



et al. 2009 1. Therefore, we also include maps of the 



image root-mean square (RMS), calculated from 
the image and coverage maps as described in Sec- 
tion [3l 

2.3. Optical Photometry and Spectroscopy 

Using the SDSS Catalog Archive Server, we 
collect the photometric information for the 97660 
sources available with data release 7 (DR7) which 
overlap with the 24 /im field of view. This data 
release covers imaging in all five optical bands, 
u',g',r',i', and z' down to r'<22.2 at a resolution 
of 1.4" in r'. After matching to our Mid-IR cata- 
log of 7603 sources, we find that 4469 of the 24/im 
galaxies have Sloan galaxy counterparts. 

We collect the Petrosian magnitudes from 
Sloan, again, replacing NGC 4921 with our own 
measurements, although in this case, the two val- 
ues are similar. 

2.3.1. Archival Spectra 

Inside the MIPS field of view, there are 346 
galaxies with SDSS spectra and redshifts between 
0.015 and 0.035. These spectra are observed with 
an aperture diameter of 3", cover the wavelength 
range from 3800 to 9200 A, and have a spectral 
resolution of 3 A in r'. Redshifts measured from 
these spectra are accurate to ^30km/s. 

We also gather spectroscopic redshifts for 869 
members from publicly available catalogs in NED 
and elsewhere in the literature. 



3 http: / / ssc.spitzer.caltech.edu /irac / calib / 
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2.3.2. WIYN Spectra 



2.3.3. Line Measurements 



To add to the Sloan spectra, we include our 
own spectroscopic observations targeting 24/im 
sources. We were awarded three nights with the 
Hydra multiobject spectrometer on WIYN in July, 
2008, however lost one night due to problems with 
telescope focusing. We observed 338 out of 529 
proposed sources which have magnitudes of r' < 
19.5, photometric redshifts of z p hot <0.2, and 
24 /im counterparts. By comparing to galaxies 
with photometric redshifts < 0.08 (see Figure [3]), 
we determine our spectroscopic sample to be 80% 
complete down to r' < 19.5, although the com- 
pleteness drops to 60% in the faintest bin (18.5 < 
r' < 19.5). The restriction in photometric redshift 
helps target a higher fraction of true spectroscopic 
cluster members, as can be seen in Figure [3j 

To obtain spectra between 3990-6795 A we in- 
tegrated for 1200 s with the blue cable and grating 
600@10.1, centered at 5380 A. This gives a disper- 
sion of 1.4 A/pixel, a spectral coverage of 2850 A, 
and a resolution of 4.6 A. 

The spectra were reduced using IDL scripts 
that trim and divide the science image by average 
dome flats, subtract the bias, perform calibration 
of the lines using arc lamps, and calibration of the 
fluxes using filtered standard stars. Cosmic rays 



were removed using the la-cosmic (van Dokkum 



2001 ) task. The scripts were run for each config- 
uration and for each night. More details on the 
reduction methods can be found in lMarleau et al.l 



(2007). 



Seventy-six of the new spectra belong to clus- 
ter member galaxies with spectroscopic redshifts 
between 0.015 and 0.035. Although many of these 
members already have published redshifts in the 
literature, with these spectra we increase the num- 
ber of member galaxies for which we can perform 
emission line diagnostics by 20%. Additionally, 
because we have targeted the MIPS 24/xm emit- 
ting galaxies in our spectroscopic survey, we in- 
crease the number of known infrared members to 
210. We gathered our spectral data in 2008, based 
on the SDSS DR6, and SDSS DR7 has since cov- 
ered many of the members we observed. Still, 
our WIYN spectra increase the number of infrared 
spectral sources by 13%. 



Optical spectra are useful threefold: First, they 
determine cluster membership; second, they allow 
for a measure of the unobscured star formation 
rate based on emission lines; and third, the emis- 
sion line ratios are a powerful diagnostic of the 
nature of line emission (either from Active Galac- 
tic Nuclei (AGN) or star formation). 

We use a custom IDL program to measure 
the emission and absorption lines for both the 
WIYN and SDSS spectra. The code makes use of 
the Markwardt package curve-fitting algorithms. 
The galaxy redshifts are determined by cross- 
correlating SDSS galaxy templates, the continuum 
is fit by hand, and the emission and absorption 
features are fit concurrently, using Gaussian func- 
tions. Close lines, such as the [Nil] -Ha triplet, are 
deblended. 

The line measurements for the 87 MIPS emis- 
sion line galaxies, for which we can access full spec- 
troscopic information, are given in Table [2] The 
full table is available in the electronic version of 
this paper, though we present the first five lines in 
both versions. 

3. Images, Spectra, and Catalogs 

We provide all of our images and photometric 
catalogs to the community online through ZR&4[[] 

The fully reduced and calibrated WIRC J, H, 
and K s images are in instrumental units. IRAC 
3.6, 4.5, 5.8, and 8.0 /im, and MIPS 24, 70, and 
160 /im images are in units of MJy sr _1 . Cover- 
age, uncertainty, weight, and RMS maps are in- 
cluded where available. The coverage maps are in 
units of number of BCDs. The uncertainty maps, 
which can be used to understand the errors, are 
computed from MOPEX accounting for Poisso- 
nian noise as well as a Gaussian component from 
the read out. The weight maps are computed in 
SWarp and are maps of the local sky background 
flux variance. RMS maps are computed on a pixel 
to pixel basis: the median value of the coverage 
map is divided by the value of the coverage map 
at each pixel, then the square root of this ratio is 
multiplied by the standard deviation of the image. 
All images discussed below are in FITS format. 



3 http:/ /irsa. ipac.caltcch.edu 
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0.01 1 1 — 1 1 — 1 

0.01 0.1 1 



spec 



Fig. 3. — Photometric redshift cut. The reliability of the photometric rcdshift estimates: spectroscopi- 
cally confirmed cluster members (shaded region) have z spec between 0.015 and 0.035 and most have z phot < 
0.08. 



Table 2 

Emission Line Measurements for Member Galaxies 



ID 


RA 


DEC 


z 


SFR JR 


Mass 


SFR Hq 


H/3 


[OIII] 


Ha 


[Nil] 




(J2000) 


(J2000) 




(M /yr) 


(10 9 M ) 


(M /yr) 




(lO^crgs^cm- 2 ) 




1 


194.757507 


26.815569 


0.0236 


0.06 


0.47 












2 


194.976395 


26.820110 


0.0268 


0.01 


3.78 


0.17 


11 ± 4 


5 ± 3 


124 ± 5 


59 ± 4 


3 


195.014603 


26.898619 


0.0196 


0.05 


31.54 












4 


194.936707 


26.954750 


0.0343 


0.33 


2.80 


0.21 


107 ± 4 


43 ± 3 


487 ± 11 


152 ± 7 


5 


194.955200 


26.974531 


0.0234 


0.26 


0.57 




116 ± 8 


77 ± 7 


374 ± 9 


81 ± 5 



WIRC - Included are intensity and weight 
maps. All images have a pixel size of 0.249". 
The image sizes are 4.3 Gbytes. 

IRAC - Intensity, coverage, and uncertainty 
maps are provided. All images are 206 Mbytes 
and have a pixel size of 1.27". 

MIPS - Intensity, coverage, and uncertainty 
maps are included. For MIPS 24/im the image 
sizes are 173 Mbytes and the pixel size is 1.27". 
For MIPS 70/im the image sizes are 18 Mbytes 
and the pixel size is 4.00". For MIPS 160/tm the 
image sizes are 2.4 Mbytes and the pixel size is 
8.00". 

We also provide the WIYN spectra we observed 
with WIYN in the region of Coma, to the NED 
spectral database. 

3.1. 24/rni Catalog 

The master catalog is based on the MIPS 24/im 
source catalog. It is constructed by finding the 
best counterpart to each MIPS source in each of 
the other bands: u', g', r', i', and z', J, H and K s 
magnitudes, and, flux densities from IRAC (3.6, 
4.5, 5.8 and 8.0 /mi) and MIPS 70/im and 160/tm. 

To discern real counterparts to the 24/im 
sources from background sources, we use the meth- 
ods of ISutherland & Saunders (19921 and Ciliegi 



et al. (2003). Briefly, this involves choosing the 
candidate with the highest reliability. The reli- 
ability is the likelihood ratio normalized over all 
possible candidates and accounts for the separa- 
tion as well as the relative magnitude distribution 
of the two data sets (24/im and counterpart). We 
consider only associations within a 5" radius of 
the MIPS 24/im source, and which have a like- 
lihood ratio greater than 0.2, see |Fadda et al.| 
(2006) for a more detailed description. The 70/im 



and 160/im images have lower resolution and we 
simply choose the candidate with the shortest pro- 
jected distance. If a 24/im source does not have 
match in a particular band, it is flagged with -999. 

Table [3] describes the machine readable 24/im 
source catalog. The first column lists the columns 
of the catalog that correspond to each entry. The 
second column lists the format of those columns, 
whether it is an integer or float, and the num- 
ber of characters. The third column gives a short 
description of each entry. Information in the ta- 
ble includes the catalog number, the J2000 RA 



and DEC position in decimal degrees, followed by 
the flux and error of the MIPS 24/im source in 
/tJy and its SNR. Next, we list the separation be- 
tween the 24/im source and its separation from 
each candidate. The SDSS u', g', r',i',z' Petrosian 
magnitudes are listed, followed by the WIRC J, H, 
and K s aperture and Petrosian magnitudes. The 
SDSS magnitudes are given in their native AB 
magnitude system. The J, H and K s magnitudes 
are in the Vega magnitude system. We list three 
aperture fluxes and the Petrosian flux for the four 
IRAC bands, and the PSF flux for MIPS 70/im 
and 160/im. We then list the photometric errors, 
respectively. This is followed by the SDSS pho- 
tometric redshift, its error, and the spectroscopic 
rcdshift. Finally several flags are listed. The first 
gives the source of the spectroscopic redshift. Val- 
ues gathered from the literature are marked with 
a T', SDSS values with a '3', WIYN values with 
'4', and NED values with '5'. The next two are 
given a value of 1 for true; the first flags close 
pairs of 24/tm sources (those within 6.1" of an- 
other source, we note that less than one percent 
have close pairs). The second flag marks the MIPS 
24/im sources associated with stars as defined by 
the SDSS star/galaxy separator. The next seven 
values list the stellarity index of, respectively, the 
IRAC 3.6 - 8.0/tm sources, followed by the WIRC 
J, H, K s sources (all calculated using SExtrac- 
tor). The final 14 flags are the close pair flags 
for the IRAC and WIRC bands followed by the 
star/galaxy flags for the IRAC and WIRC bands. 

4. Results 

Our aim is to derive the specific star formation 
rate (sSFR) for Coma galaxies by taking advan- 
tage of the large baseline of multi- wavelength data 
available. As a first step, we classify the galaxy ac- 
tivity as either star formation, or as contaminated 
by AGN. We employ several diagnostics, including 
optical emission line ratios, IRAC colors, and the 
Radio Far-Infrared (radio-FIR) correlation. Next, 
we build the SEDs based on up to 15 photometric 
datapoints and match these to empirical template 
galaxy spectra from which we measure the total 
infrared luminosity, and galaxy stellar mass. We 
find that galaxies with specific SFRs > 0.1 Gyr -1 
are blue, and that the majority are dwarfs, having 
mass<lO 9 M . 
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Table 3 

MIPS 24 fiM Source Catalog Columns 



Column 


Format 


Description 


1 


i6 


Catalog Number 


2 


f'14.6 


RA (J2000) 


3 


f'14.6 


DEC (J2000) 


4 


fl0.3 


MIPS24 (^ijy) 


5 


fl0.3 


MIPS24 error (/J.Jy) 


6 


f9.3 


SNR 


7 


f9.3 


MIPS24-SDSS Sep (") 
MIPS24-IRAC Sep( // ) 
MIPS24-NIR Sep(") 
MIPS24-70 Sep(") 
MIPS24-160 Sep(") 
MIPS24-Z Sep(") 


8 


4xf9.3 


12 


3xf9.3 


15 


f9.3 


16 


f9.3 


17 


f9.3 


18 


fl0.3 


u (mag) 


19 


fl0.3 


g' (mag) 


20 


fl0.3 


r' (mag) 


21 


fl0.3 


i' (mag) 
z' (mag) 


22 


f'10.3 


23 


2xfl0.3 


J ap 3.5" , Pctrosian (mag) 


25 


2xfl0.3 


H ap 3. 5", Pctrosian (mag) 


27 


2xfl0.3 


K s ap 3.5" ,Pctrosian (mag) 


29 


4xfl0.2 


IRAC1 ap 4 ,ap 6 ,ap 12' .Pctrosian (/iJy) 


33 


4xfl0.2 


IRAC2 ap 4",ap 6",ap 12" ,Petrosian (/iJy) 


37 


4xfl0.2 


IRAC3 ap 4",ap 6",ap 12" , Pctrosian (pjy) 


41 


4xfl0.2 


IRAC4 ap 4",ap 6",ap 12" , Pctrosian (pjy) 




flO 2 


IvllroiV {fJ/jy ) 


46 


fl0.2 


MIPS160 (/jjy) 


47 


f'10.3 


u error (mag) 


48 


f'10.3 


g' error (mag) 


49 


f'10.3 


r' error (mag) 


50 


fl0.3 


i' error (mag) 


51 


f'10.3 


z' error (mag) 

J error ap 3.5" , Pctrosian (mag) 


52 


2xfl0.3 


54 


2xf'10.3 


H error ap 3.5" , Pctrosian (mag) 


56 


2xfl0.3 


K s error ap 3.5" , Pctrosian (mag) 

IRAC1,IRAC1 error ap 4",ap 6",ap 12" .Petrosian (p.ly) 


58 


4xf'10.2 


62 


4xfl0.2 


IRAC2,IRAC2 error ap 4",ap 6",ap 12" , Pctrosian (jiJy) 


66 


4xfl0.2 


IRAC3,IRAC3 error ap 4",ap 6",ap 12" ,Petrosian (fiJy) 


70 


4xfl0.2 


IRAC4,IRAC4 error ap 4",ap 6",ap 12" , Pctrosian (>Jy) 


74 


f'10.2 


MIPS70 error (/jjy) 


75 


f'10.2 


MIPS160 error(/iJy) 


76 


f'10.2 


SDSS Photz 


77 


f'10.2 


SDSS Photz error 


78 


f'10.2 


Spec z 


79 


f'10.2 


Spec z source l:LIT,3:SDSS,4:WiYN,5:NED 


80 


f9.3 


MIPS24-Z Sep(") 


81 


i2 


MIPS24 separation flag 


82 


i2 


MIPS24 star/galaxy flag 


83 


7xf7.3 


IRAC1-4,J,H,K stellarity index 


90 


7xi2 


IRAC1-4,J,H,K close pair flag 


97 


7xi2 


IRAC1-4,J,H,K star/galaxy flag 
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4.1. AGN Contamination 

AGN activity can mimic many of the same 
emission characteristics as young starbursts, in- 
cluding bright MIPS detections. Using several ad- 
ditional diagnostics, we show that the contamina- 
tion from AGN is low. 

4- 1.1. Optical Emission Lines 

Of the 210 known cluster members with MIPS 
24/im emission, 125 have spectra. Of these last, 
38 are passive galaxies with no measurable emis- 
sion lines. Table [2] presents the emission line mea- 
surements of the remaining 87 used to construct 
the BPT diagram ( |Baldwin et al.||1981[ ) of Fig- 
ure [2] The dashed lines separate the star forming 
and AGN regions oflOsterbrock & Ferlandl (I2006I) , 



and the curved lines show separations between star 



forming galaxies from the SDSS (Kauffmann et al. 



2003), composite sources, and a theoretical up- 
per limit for sources that can be described by star 
forming models (Kewley & Dopita 2002). Only 
four are Seyfert galaxies, and 14 LINER. Most 
emission line galaxies fall in the star forming re- 
gion, or in the transition region where star forma- 
tion is still a dominant process. 

4.1.2. IRAC Colors 

Stars, galaxies, and broad line AGN can be dif- 



ferentiated using their IRAC colors (Eisenhardt 



et al.||2004| |Lacy et al.|[2004l |Sajina et al.||2005 



Donley et al. 2008). AGN are well separated be- 
cause the power law of AGN emission is more red 
than the galaxy spectrum in the 3.6/xm channel, 
and because there is a lack of PAHs in AGN. |Lacy| 
etal] ([20041) plot the IRAC colors as (8.0/4.5) vs 



(5.8/3.6) which separates objects with blue con- 
tinua from objects with red continua; AGN tend 
to be red in both IRAC color filters. 

We plot the IRAC colors of our MIPS selected 
galaxies in Figure [5j The majority of bright IR 
sources avoid the AGN region, although four may 
be classified as AGN. The predicted IRAC colors 
from template galaxies used in our SED fitting, at 
the redshift interval of the Coma Cluster, are also 
shown in the figure. As can be seen, the template 
AGN (triangles) occupy approximately the same 



4-1.3. Radio Far-Infrared Correlation 

Using the radio-FIR correlation for star forming 



galaxies of Condon ( 1992 ) , we can determine the 



presence of radio-excess AGN. We use the cata- 
log of 1.4GHz radio fluxes for Coma galaxies from 



Miller et al. (2009) and match the radio galax- 
ies to our MIPS 24/zm and 70/zm sources. The 
1.4GHz resolution is similar to that of the 24/im 
resolution, and we match our sources purely on 
proximity in RA, Dec, and redshift. As shown 
in Figure |6j sources well below the relationship 
harbor a radio excess, this is attributed to AGN 
emission. From Appleton et al.| ( 2004 1 , the limit- 
ing ratios are q 2 4= log (S 2 4 /1 m/Si.4Gffz)~ 1 ± 0-27 
and q 70 = log (S 70 „ ro /Si.4GH*) ~ 2.30 ±0.16. We 
include a small color correction of -0.017 for 24/im, 
and the 70/im relation includes the updated aper- 
ture correction of +0.15 which were not applied 



et al. 


2004 


Frayer et al. 


2006) 



the sources do not have this radio excess, and are 
therefore dominated by star formation. 

4- 1-4- AGN Contamination is Low 

Using the three methods above, as well as the 
fits to SEDs discussed in Section [4.2.11 we identi- 
fied only 30 possible AGN. One source was identi- 
fied from three of the diagnostics: the Radio FIR 
correlation, the IRAC colors, and the SED fit; 
three sources were identified by 2 methods (Radio 
FIR+BPT:1, BPT+IRAC:1, IRAC+SED:1); but 
most sources were identified by only one method 
(Radio FIR: 7, BPT: 14, IRAC: 3, SED:2). There- 
fore, star formation is the most important contrib- 
utor to galaxy activity in the Coma galaxies of our 
sample. This represents a conservatively large es- 
timate of the contribution of AGN activity, as the 
bolomctric luminosity in these systems is not nec- 
essarily dominated by the AGN. For example, as 
discussed in Section |4~2"l the best fit SEDs immedi- 
ately identify dominant AGN - only 4 of 210 con- 
firmed Mid-IR members. Furthermore, |Veilleux] 



locus as defined by Lacy et al. (2004) 



et al. (2009) derive a diagnostic using the ratio 
of 25//m/60//m flux for ultraluminous IR galaxies 
and Palomar-Green QSOs (PGQSOs). Using this 
diagram, only eight of our possible AGN corre- 
spond to > 50% AGN contribution, and 22 < 50% 
AGN contribution. 
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Fig. 4. — BPT diagram. Ratios of nebular emission lines can be used to separate optical AGN from clear 
star forming galaxies. The dashed lines separate the star forming and AGN regions of Osterbrock fc Ferland| 
( |2006[ ) and the curved lines separate star forming, composite, and AGN sources (Kauffmann et al. 2003). 
Plotted, are the cluster members which contain measurable line emission. The majority fall in the starburst 
and transition regions of the diagnostic diagram, with few clear AGN. Larger dots represent higher specific 
star formation rates, based on SED fitting, and assuming all the luminosity comes from star formation (see 
the discussion in Section [4.2.31 



o 

CO 

M 
O 




0.5 

Log(5.8/3.6) 



Fig. 5. — IRAC colors. Stars, ellipticals, and spirals can be separated using IRAC colors. We plot MIPS 
selected members which have detections in all four IRAC bands, following Lacy et al. ( |2004 ). Black dots 



are for passive galaxies, small filled dots are low and star forming galaxies, and the larger filled dots are 
starbursts. Open shapes represent the IRAC colors from several templates: ellipticals (diamonds), spirals 
(circles), AGN (triangles), starbursts (squares). Very few of the spectroscopically confirmed members are 
inside of the dashed polygon, the AGN region of the plot. 
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Fig. 6. — The radio-FIR correlation. The plots show the logarithm of the observed flux density ratio 
(q) versus the infrared flux density for all spectroscopic members with both a MIPS and radio detection. 
Sources with MIPS 24/xm emission are on the top, and with MIPS 70 /xm emission are on the bottom. The 
solid grey bar indicates the best-fit radio-FIR correlation and its 3a error. The grey hashed area under the 
correlation band is expected to be populated only by radio-excess AGN. 
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4.2. Obscured Star Formation 

At this point, we have used spectra to iden- 
tify cluster members, and are confident that there 
is little contamination from AGN. We now quan- 
tify the amount of star formation by deriving the 
specific star formation rates based on the total in- 
frared luminosity and the galaxy stellar mass. The 
total infrared luminosity is produced from dust 
heated by hot young stars, and as such provides a 
measure of obscured star formation. 

4-2.1. Spectral Energy Distributions 

We build the measured SED for the 210 con- 
firmed 24/im emitting galaxies, using up to 15 
photometric bands. These SEDs are compared 
to semi-empirical templates from |Polletta et al.| 
( 2007 ) , which are particularly well adapted to fit- 



ting infrared emitting galaxies. There are several 
templates for Elliptical, Spiral, Starburst galaxies, 
and AGN. These are based on results from the 



GRASIL (Silva et al. 1998) code combined with 



photometric data from IR galaxies and AGN. The 
best-fit is the template with the lowest chi-square. 
Figure [7] shows a few examples of confirmed cluster 
members which have photometry in the represen- 
tative bandwidths (optical, Near-IR and Mid-IR) 
and redshifts between 0.015-0.035. The total in- 
frared luminosity is found by integrating the flux 
density between 8-1000/iin of the best fit template. 
Only 4 of the 210 cluster members correspond to 
a best-fit template which is completely dominated 
by an AGN. 

4-2.2. Galaxy Stellar Mass 

To determine the galaxy stellar mass, we use 



the stellar model templates of Maraston (2005) 



with a Kroupa initial mass function (Kroupa 



2001), adopting a solar metallicity. These models 



describe the stellar emission from the UV to the 
Near-IR and we limit the fits to the optical and 
Near-IR portion of our observed SED, up to only 
IRAC 3.6 and 4.5/xm. We modify the template 



SEDs in accordance with the;Calzetti et al. (2000) 



extinction law, allowing E(B-V) to vary between 
and 1. We then fit our observed SEDs to the 



host of available Maraston ( 2005 ) templates, and 



retrieve the stellar mass from the best-fit model. 
Photometric points from our Near-IR photom- 
etry, and those from the IRAC bands help to 



tightly constrain the SED fit around the stellar 
mass bump. 

Table [2] gives the calculated stellar mass for the 
210 MIPS members. 

4-2.3. Specific Star Formation Rates 

We convert the total infrared luminosity to a 
star formation rate (listed in Tab le [2| using the 
Kennicutt relationship ( Kennicuttp 998). 

The specific star formation rate is simply the 
SFR divided by the stellar mass. It is a useful 
parameterization since Coma is host to a variety 
of galaxies. For example, a dwarf starburst galaxy 
may have a much smaller SFR than a large passive 
elliptical, yet a higher value of sSFR. In most of 
our discussion, we present the ratio f s &, defined as 
the sSFRxt, the timescale over which the galaxy 
is assumed to form stars at the currently observed 
rate. We take r = 100 Myr following |Fadda et al 
(2006). Starburst galaxies here are those with 



t* >o.i. 

Figure [8] shows the specific star formation rates 
with the galaxy flux at each of the Mid-IR wave- 
lengths. A range of 24/im luminosities give high 
specific star formation rates, suggesting that fit- 
ting the total SED, which finds both the mass, 
as well as the star formation rate, may provide a 
break in the degeneracy that exists when deriving 
star formation rates purely from the 24/um lumi- 
nosity. 

We separate Figure [8] into four regions of sSFR. 
The passive galaxies have sSFR < 5xl0~ 4 Gyr -1 , 
low star formers have 5xl0 -4 < sSFR < 5xl0 -3 
Gyr -1 , normal star forming galaxies have 5xl0 -3 
< sSFR < 0.1 Gyr -1 , and starbursts have sSFR > 
0.1 Gyr -1 . Many of the MlPS-detected galaxies 
are passive, and the IR emission is likely a result 
of dust expelled by the large number of stars at 



later stages of their evolution (Temi et al. 
Calzetti et al]|2010[ ). 



2008 



The best-fit SED galaxy type, as given by the 
empirical SED, is plotted as a histogram with spe- 
cific star formation in Figure [8] The Elliptical, SO, 
and Spiral galaxy types correspond well to the ex- 
pected behavior of the designations we set based 
on the sSFR, with passive (< 5xl0 -4 Gyr -1 ) cor- 
responding to ellipticals. The few AGN are also 
displayed. 
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Fig. 7. — Sample best-fit template SEDs. Measured photometric points are shown on top of the best-fit 
Polletta (2007) semi-empirical template SED. 
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Fig. 8. — Specific star formation rates. The Mid-IR luminosity of cluster members is plotted as a 
function of the specific star formation rate, derived from SED fitting. There are four regimes: Passive, 
Low star formation, Normal star formation, and Starburst. Grey points indicate member galaxies within 
0.3 degrees of the cluster X-ray centroid. Top - The 24/xm luminosity. Second - The 70/im luminosity. Third 
- The 160/im luminosity. Bottom - The fraction of galaxy types as a function of specific star formation rate. 
The galaxy types are those of the best fit templates from Polletta (2007). Galaxies are classified as ellipticals 
(thin red line), SO (thick black line), spirals (thick blue dashed line), and AGN (thin green dashed line). 
Ellipticals dominate at low sSFR where red galaxies are the majority (see Figure |9|. Only four sources have 
best-fit templates corresponding to AGN. 
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5. Discussion 

We now turn to a discussion of the star forma- 
tion rates based on SED fitting. In examining the 
galaxy colors, we find that the most red galaxies 
have the lowest sSFRs, and that a few galaxies 
with normal specific SFRs are also red. 

We compare the sSFRs based on SED fitting, to 
values based purely on MIPS 2Afim observations, 
confirming earlier results and extending to lower 
stellar mass. Additionally we find that fits in- 
cluding 70/im data estimate slightly higher sSFRs 
than those from MIPS 24^m observations alone. 

We compare the obscured sSFRs to extinction- 
corrected star formation rates derived from Ha 
emission lines. As we will see, the extinction cor- 
rected unobscured values are able to completely 
recover the obscured values. 

Finally, we describe the sSFR as a function of 
local galaxy density, finding a higher fraction of 
star forming galaxies at the cluster outskirts, than 
within the cluster core. 

5.1. The Color and Mass of Star Forming 
Galaxies 

Figure [9] shows the galaxy optical color and 
stellar mass, as a function of sSFR. The star- 
burst galaxies are blue, with (g'-r')< 0.7, and are 
dwarfs, with mass < 10 9 M Q . Most sources with 
a specific star formation rate less than 0.005 are 
red and have a large stellar mass. These are the 
red-sequence galaxies that dominate the cluster 
core and the FIR emission is likely caused by a 
large population of old red giant branch stars (also 
known to emit at MI PS 24/zm; |Temi et al.] ( |2008| ) ; 
Calzetti et al] (|2010[)). 



5.1.1. No Evidence for Dusty Starbursting Dwarfs 

One might expect that if a large population 
of dusty dwarfs existed, some low mass galaxies 
should be red by virtue of their dust extinction, 
rather than because of age. There are a few red 
dwarfs with normal star formation, but we do not 
detect any starbursting red dwarfs. 

Because dust-obscured starburst galaxies should 
be optically fainter than their unobscured, blue, 
counterparts, it may be that our spectroscopic 
data is too shallow or incomplete to detect these 
fainter dwarfs. To investigate this issue, we 



present Figure 10 which makes use of a deeper, 
complete catalog of sources based on photomet- 
ric redshifts. The figure plots the color magnitude 
diagram of MIPS detected members. The spectro- 
scopic members (filled black circles) are compared 
to the photometric members (open circles; defined 
as those with z p / lot <0.08). If it was true that 
our spectroscopic limit introduced a bias against 
dusty red star forming dwarfs, than the photomet- 
ric members, which are much fainter, would show 
a sample of low-mass red galaxies. But, as the fig- 
ure shows, there does not appear to be a missing 
population of red 24^m sources just beyond our 
spectroscopic limit of r'^19.5. There is a caveat. 
Unfortunately, it is not possible to interpret the 
data below r'~20.5 as the photometric redshifts 
become unreliable. Nevertheless, the lack of dusty 
red dwarfs does not appear to be caused by the 
spectroscopic limit of our study. 

5.2. Comparison of Star Formation Rates 
from other Methods 



5.2.1. Obscured 



Unobscured Rates 



To derive the star formation rate for the un- 
obscured component we measure the Ha emission 



line and use the conversion of Kennicutt ( 1998 ) 



The Ha emission was observed through fixed- 
diameter fiber spectrographs, but Coma galaxies 
are relatively nearby, so many of the galaxies are 
larger than the ^3" fiber diameter. Thus, we ap- 
ply an aperture correction. Using SExtractor on 
the SDSS r' image, we retrieve the isophotal area 
for each source. We assume that the star forma- 
tion is spread throughout the galaxy area, thus 
the ratio of the galaxy area to fiber area gives the 
aperture correction. Corrections for the largest 
Coma members can be substantial, up to a factor 
of ten, however, the median correction is ^2. 

Star formation rates based on optical emission 
lines suffer greatly from extinction by surrounding 
dust, and so a second important correction we ap- 
ply to the Ha flux measurement is an extinction 
correction. In most sources, we can measure the 
H/3 emission line which we use to find the Balmer 
decrement and calculate the amount of extinction. 
We calculate a significant E(B-V) with a median 
of 0.55 and standard deviation of 0.42. Extinction 
and aperture corrected star formation rates from 
the optical emission lines are given in Table [2j 
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Fig. 9. — Galaxy Color and Mass. The color and mass of Coma members are plotted as a function 
of the specific star formation rate as measured by SED fitting. There are four regimes: Passive, Low star 
formation, Normal star formation, and Starburst. Top - The optical (g'-r') color. The grey points highlight 
galaxies within 0.3 degrees of the cluster core. - Bottom The best-fit SED mass is shown, with a grey line 
showing our separation of dwarf galaxies with Mass < 10 9 M Q . Red-sequence galaxies (see Figure 10) arc 
outlined in red. All of the starburst galaxies have blue colors with (g'-r')< 0.7, and are dwarfs. 
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Fig. 10. — Optical color magnitude diagram. The color magnitude diagram of Coma members with 
MIPS detections. Photometric members (with z p hot< 0.08) are plotted as open circles and show the depth of 
the photometry. Spectroscopic members are filled black points. Large dots represent the dwarf galaxies. Our 
spectroscopic limit of r'=19.5 is delineated with a vertical dotted line. There are no low mass red galaxies in 
the spectroscopic sample. There is also no evidence from the photometric sample that a significant population 
of red members exists just below the spectroscopic threshold. Therefore, the lack of dusty red dwarf galaxies 
is not because of our spectroscopic limit or incompleteness. 
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SED SFRs Reproduce Extinction Corrected Ha 
Values 

We can test how well the star formation rates 
based on total IR luminosities compare to values 
measured from corrected Ha emission line fluxes. 
In Figure [TT] we show the ratio of the SFR de- 
rived from the optical emission line flux to that 
from SED fitting. 

The aperture correction is more important for 
the larger galaxies, and once performed brings a 
uniform ratio of unobscured to obscured star for- 
mation for galaxies with significant SFRg^ (Fig- 
central panel). 



ure 



11 



Including the extinction 
correction, performed on a source by source ba- 
sis, the two rates approach each other. The me- 
dian value of LogpFRfla/SFRss^)] is 0.0, and the 
standard deviation is 0.5. This is an important 
point, for the sample of galaxies that have both 
emission lines and MIPS detections, the full star 
formation rates of the optical galaxies can be re- 
covered after correcting for the extinction using 
optical emission lines. We suggest that this is be- 
cause the same dust is causing extinction of the 
optical lines, as is producing the Mid-IR emission. 

5.2.2. SFR Estimates from 24/J-m Flux only 

Many studies have used the 24/im emission 
alone as a proxy for total infrared luminosity 
([Rieke et al.||2009||Bai et al.|2006||Mahajan et~aT 
2010 ). Certainly the SED fitting technique is more 



robust, as it is based on up to 15 photometric 
points. That the method is more accurate, is 
still open for debate, however, recent results from 
the Herschel space telescope support this. These 
new data show how star formation rates based on 
24/im data underestimate rates derived by mea- 
suring longer wavelength fluxes, in 40% of the 



cases ( |Rawle et al. 2010). Our data also show 
a region where the 24/im data underestimates the 
rate calculated from the full SED fitting. In Fig- 
ure 12 our results agree with the relationship of 



Rieke et al. (2009) (the dashed line), out to the 



high flux end. Beyond 4mJy, the 24/im under- 
estimates the IR luminosity as measured by the 
best-fit SED. 

5.2.3. Specific SFR from (24fJ,m - 2!) Color 
In this section, we show that our results con- 



who have recently explored the specific star for- 
mation in Coma cluster galaxies. These authors 
have used the (24/im - z') color as a proxy for 
quantifying sSFR. Our method relies on selecting 
galaxies based on the physical criteria of stellar 
mass and SFR. 

The first panel of Figure [13] shows that the 
usage of (24/im - z') recovers normal star form- 
ing and starbursts, galaxies with specific SFR > 
0.01 Gyr -1 , although, their color cut of -6 also in- 
cludes contamination by many lower sSFR galax- 
ies. There is also a population of galaxies with low 
and normal specific star formation rates (between 
0.0005-0.1 Gyr" 1 ) which lie above the (24/im - z') 
color cut of -6. 

One of the results of our work is that the in- 
frared star forming and starburst galaxies are blue, 
having (g' - r') < 0.7. The central panel of Fig- 
ure 13 confirms that the (24/im - z') > -6 selects 
the same blue and red populations based on the 
optical color cut. 

The right panel of Figure [13] shows that we 
are able to probe deep into the low-mass regime, 
measuring galaxies with stellar masses as low as 
10 7 5 M M . The vertical dotted line shows the limit 



firm and extend those of Mahajan et al. (2010), 



of the Mahajan et al. (2010) paper at a mass of 
10 8 - 5 M Q . As we saw in Figure 9j the low mass 
dwarf galaxies are important as the highest spe- 
cific sSFR galaxies, the starbursts, all have masses 
lower than 10 9 M Q . 

5.3. Environment 

Most of the imaging extends beyond the virial 
radius of Coma, parametrized by r 2 oo, the char- 
acteristic radius for the cluster. We can therefore 
compare the galaxies in the cluster core, 0.2xr 20 o, 
to those on the edge of the cluster center, beyond 
half of r 2 oo- Fo r Coma, r 2 oo has b een measured to 
be l.Sh^Mpc ( |Geller et al|l999| ). At the cluster 
rcdshift, this translates to r 2 oo= 88'. 

5.3.1. Location 

Figure [14] shows the positions of the spectro- 
scopically confirmed 24/im emitting cluster mem- 
bers on the sky. Most dwarf galaxies are starbursts 
or normal star forming galaxies, and few exist in 
the cluster core. 

The suppression of star formation in the clus- 
ter core is also is consistent with lCastander et al.l 



21 



No Corr 





Ap Corr 





Ext And Ap Corr 



wL 



wL 



0.001 0.01 0.1 



10 0.001 0.01 0.1 



10 0.001 0.01 0.1 



10 



SFR 5ED [M s yr- 



SFR SED [M e yr-'] 



SFR SED [M yr-'] 



Fig. 11. — Unobscured versus obscured star formation rates. The ratio of unobscured to obscured 
star formation rates as a function of obscured star formation rate. Left No correction is applied. Center 
- Only the aperture correction is applied. Right - Both aperture and extinction corrections are applied. 
Much of the optical emission in these nearby Coma galaxies is beyond the fiber diameter, so the aperture 
corrections are important (central panel). Extinction corrections raise the mean value of the sample. Once 
the extinction correction is applied, the Ha based SFR matches the rate based on the total IR luminosity 
of the best-fit template. Any source that has been flagged as possible AGN in Section 4.1 is shown as a 
triangle. There are few, and they inhabit no particular region of this diagram. 
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Fig. 12. — Total infrared luminosity versus 24/im flux. The total infrared luminosity from SED fitting 
is plotted against the measured 24/im flux. The dashed line shows the relationship developed by |Rieke et al.| 
(2009). Sources that have been flagged as possible AGN are shown as triangles. The overall agreement is 



good, however, beyond 4 mJy the 24/im flux relationship underestimates the IR luminosity as measured by 
the best-fit SED. 
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Fig. 13. — The (24/im - z') color. Left - We compare the specific star formation rate of cluster galaxies 
using two parameterizations. On the y-axis we plot the (24/im - z') color, and on the x-axis, we plot the 
sSFR as determined from our SED-fitting technique. The vertical dotted lines define Passive, Low, Normal 
Star forming, and Starburst galaxies as in the text. The larger red points highlight red-sequence galaxies. 
Center - (24/im - z') color versus optical color, plotted to a (g' - r') of 1.0. Here, and in the right panel, 
dotted lines show limits from Mahajan et al.| ( |2010 ), the horizontal dashed line is at (24/im - z') — -6, the 
separation used in Mahajan et al. ( |2010 ), and larger black points indicate galaxies in the low-sSFR regime 
according to our SED fitting. Right - (24/im - z') color versus galaxy stellar mass. We confirm the results of 
the Mahajan et al. (2010) study, and extend to lower stellar mass. This last is enabled by our deep Ncar-IR 
photometric points. 



(2001) who found a dearth of active galaxies in 



the core of the Coma cluster based on SDSS spec- 
troscopy. 

5.3.2. Local Galaxy Density 

We compute the galaxy density using the posi- 
tions of 1557 red-sequence galaxies from the SDSS 
with 12<r'<20, using a 2D wavelet transform as in 



Fadda et al. (1998). The red-sequence is the same 
as shown in Figure 10, where (g'-r') = 0.036 xr' - 



1.320, and we include a range of ±0.07 in (g'-r') 
color. The units of galaxy density are thus in num- 
ber of red-sequence galaxies (Nrs) per square de- 
gree, represented as contours in Figure 14 The 



core and well-studied group NGC 4839 at (194.35, 
27.5) are easily detected. A previously undetected 
group at (194.1, 27.7deg) is also detected along the 
direction of the filament toward Abell 1367. 



The top panel of Figure 15 reveals a dearth of 
starburst galaxies in the cluster core, where the 
galaxy density n^js<300. In this figure, we aim 
to examine the star forming galaxies with both 
high, and low sSFRs. Thus, in the histograms 



we have excluded passive galaxies, those on the 
red sequence and having masses >10 10 M Q . We 
also removed the 4 galaxies with SEDs best fit by 
AGN. In the second panel, the number of galaxies 
with a high sSFR (blue solid line) or a low sSFR 
(red-orange dashed line) is normalized by area. N 
increases toward the cluster core for both popula- 
tions. The number of star forming galaxies with 
high specific SFR (>0.005 Gyr -1 ) is lower than 
the number of low sSFR galaxies with the most 
differentiation in the cluster core. Even when nor- 
malized by the number of red sequence galaxies 
(third panel), the modestly star forming galaxies 
dominate over normal and SB galaxies, for all den- 
sity bins. 

Only a minority of passive galaxies are detected 
in the infrared, therefore, the number of red se- 
quence galaxies (which may or may not have sig- 
nificant infrared emission) , is much larger than the 
number of galaxies which we have defined as be- 
ing passive from their low specific star formation 
rates. The number of high specific SFR galax- 
ies normalized by the number of low specific SFR 
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Fig. 14. — Spatial distribution of star formation. The position of each spectroscopically confirmed 
MlPS-emitting members on the sky is shown. Red points are passive galaxies, black are low specific star 
formers, blue are normal star forming and cyan are starburst. The dwarf galaxies are shown as large points, 
and galaxies with masses greater than 1O 9 M are smaller points. Seyfert galaxies are shown as black triangles. 
The grey circle has a radius of 0.3 degrees and highlights the core region of the cluster, the black rectangle 
marks the MIPS field of view. Contours show the galaxy density, i.e., the number of red-sequence galaxies 
per square degree. As can be seen, the core has a dearth of starburst galaxies. 
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Fig. 15. — Local galaxy density. Top - The specific star formation rate for cluster members as a function 
of galaxy density. The galaxy density is in units of number of red-sequence galaxies per square degree. A 
density of 900 corresponds to the core of Coma and is indicated with a dotted line. The filled blue circles 
represent galaxies with normal SF and as well as those currently undergoing a SB. The open red-orange 
circles represent galaxies with low and no star formation. The largest circle size showcases the Dwarf galaxy 
population, and the smallest circles are the massive rcd-scqucncc galaxies (as described in the text). For 
these last, we assume the MIPS emission is the old stellar population, and therefore, they are omitted from 
the histograms below. Second - The number of galaxies per degree squared in bins of local galaxy density. 
The red-orange dashed line corresponds to low SF and passive galaxies, the blue solid line corresponds to 
normal and SB galaxies. Third - The number of galaxies normalized by the total density in the bin. There 
are fewer active galaxies throughout, with the largest fractional difference being in the core of the cluster, 
where passive galaxies dominate. Bottom - The fraction of galaxies with high (sSFR>0.005 Gyr -1 ) to low 
(sSFR<0.005 Gyr -1 ) specific SFRs. The fraction of high-to-low star forming galaxies decreases toward the 
Coma core. We show Poissonian errors on the histograms. 
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galaxies is shown in the bottom panel. Indeed 
the trend is the same as for the solid blue line of 
Nss/N^s, with a progressive decrease toward the 
cluster core, but the numbers are quite different. 
We detect only ~10% of the red sequence galax- 
ies as infrared, low star forming galaxies outside 
the cluster core, and roughly ^30% in the high- 
est density regions of the core. This is consistent 
with a larger fraction of red sequence galaxies hav- 
ing been recently quenched in the cluster core. As 
a caveat, we reiterate that the observations of the 
Coma do not extend much past the virial radius, 
so we are not sampling the filament that connects 
to Abell 1367. 

The Distribution of Specific SFRs 

To further explore the effect of galaxy density, we 
examine the distribution of specific SFRs. The 
left hand side of Figure [16] shows the distribution 
of specific star formation rates for the core (dashed 
line; density >900), intermediate regions (dotted 
line; 300< density <900), and the low-density out- 
skirts (solid line; density <300). We again exclude 
AGN and large red galaxies, so the number of core 
galaxies is small. The right hand side of the fig- 
ure shows the same, but normalized by the total 
density in each bin. The core galaxies have a dis- 
tribution of sSFR skewed to lower values, peaking 
near 5xl0 -4 Gyr -1 , which, according to our SED 
fitting and Figure|8j corresponds to where the frac- 
tion of elliptical galaxies begins to decrease. In 
other words, the core is dominated by a larger 
number of galaxies with low sSFR, as expected. 
The low density regions have a Gaussian distri- 
bution around 5xl0 -3 Gyr -1 sSFR characteristic 
of quiescent (SO) and normal star forming galax- 
ies. The intermediate density regions show both 
these populations, this may be a zone of transition 
between the two types. 

6. Summary and Conclusions 

We have exploited a plethora of archival as well 
as original photometric and spectroscopic observa- 
tions in order to study the specific star formation 
of infrared member galaxies of the Coma cluster. 
Of particular importance for SED fitting is the 
new Near-IR photometry we collected at Palomar, 
as well as the inclusion of the longer wavelength 
MIPS observations. We matched the 24/im photo- 



metric catalog to spectroscopy gathered from the 
SDSS and to our own Hydra/WIYN observations 
to confirm 210 cluster members. We show that 
the contribution of AGN in the sample is low, and 
build spectral energy distributions which we fit to 
templates in order to derive the total infrared lu- 
minosity and galaxy stellar mass. 

We compare these estimates of obscured star 
formation to unobscured star formation rates de- 
rived from optical emission lines, finding that the 
two values agree once the optical spectra have 
been corrected for aperture size and extinction. 
This provides confidence to studies restricted to 
ground based observational methods only. For the 
filament and outskirt regions, where background 
contamination is high, the Mid-IR remains a good 
tool for initially identifying star forming galaxies. 

By further exploring the optical colors and spe- 
cific star formation rates based on total infrared 
luminosities, we find that the starburst galaxies 
are blue dwarfs. This confirms the results of IMa-l 



hajan et al. (2010) and extends their results to 



lower mass galaxies using a more robust probe of 
galaxy SFR and a more accurate measure of mass. 

Additionally, we find that the star formation 
rates based on MIPS 24/im data alone represent 
well the values derived from fits to the total SEDs, 
except at high fluxes, where rates are slightly 
underestimated compared to those based on the 
SED fitting technique. We conclude that multi- 
wavelength photometry, in combination with fit- 
ting to template SEDs provides the most robust 
estimate of specific star formation rates. This is 
because the infrared hump is better constrained 
by values beyond the maximum. It is also a good 
galaxy identifier, as AGN are immediately spotted 
from their SEDs. 

Confirming much of the recent work on Coma 
and nearby clusters in general, we also find that 
most of the star forming galaxies lie outside of the 
core. In particular, we have found that the ra- 
tio of high sSFR galaxies to low sSFR galaxies is 
lowest at the cluster core. We have also found 
that the distribution of the sSFR is skewed to- 
wards lower rates for the galaxies in high density 
regions, even after removing luminous red galax- 
ies. This proves the core to be dominated by a 
larger number of low-SF galaxies. Furthermore, 
the intermediate density regions show a bimodal 
distribution in sSFR, with peaks at the same rates 
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Fig. 16. — Distribution of sSFR. The distribution of specific star formation rate is shown for the three 
density regions considered in Figure [15] The core, with density >900 is the dashed line, the intermediate 
regions are shown with a dotted line, and the lowest density regions with <300 are shown as the solid line. 
Left - The histogram of sSFR. Right - The distribution of sSFR for the three density regions, normalized by 
the total Nrs in each bin. Large red galaxies, which likely have MIPS emission only from a population of 
old stars, are excluded from the distributions as are the 4 AGN identified from SED fitting. The galaxies in 
the high density core show a distribution of specific SFRs skewed to lower values. 



of both the high and low density regions. We sug- 
gest therefore that these intermediate density re- 
gions are the transition zone between the low and 
high sSFR galaxies, in line with the popular no- 
tion that star formation is quenched in galaxies 
before they enter the cluster core. 

Finally, we provide to the community fully re- 
duced images in J, H, K s , IRAC 3.6, 4.5, 5.8, and 
8.0/im, and MIPS 24, 70, and 160^m through 
IRSA and 338 new spectra of 24^m sources 
through the NASA Extragalactic Database. 
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